


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1970-09 


An analysis of the sets of rank distribution in 
an hierarchical organization 


Armacost, Robert Leo 


Monterey, California ; Naval Postgraduate School 


http://ndl.handle.net/10945/14880 


Downloaded from NPS Archive: Calhoun 


: Calhoun is the Naval Postgraduate School's public access digital repository for 
/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 








United States 
Naval Postgraduate School 





THESIS 


AN ANALYSIS OF THE SETS 
OF RANK DISTRIBUTION 
IN AN HIERARCHICAL ORGANIZATION 


by 


Robert Leo Armacost 





September 1970 


This document has been approved for public re- 
Lease and sake; its distribution 416 unloncted. 


LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIF. 93940 


LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTERLY, CALIF. 93940 


An Analysis of the Sets of Rank Distribution 


In an Hierarchical Organization 
by 


Robert Leo Armacost 
Lieutenant, United States Coast Guard 
B.S., United States Coast Guard Academy, 1964 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN OPERATIONS RESEARCH 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 1970 





ABSTRACT 

A Markov type model for studying rank distribution in 
an hierarchical organization is examined. Various sets of 
rank distributions are defined and their properties discussed. 
A computer aided test is developed for testing a given rank 
distribution to determine if it is an element of a particular 
set. For a distribution not in the Steady State set, a test 
is developed to determine if the distribution can be returned 


to in m transitions or steps. 
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I. INTRODUCTION 

Several authors have studied models for analyzing rank 
distributions in an hierarchical organization. Bartholomew 
[1967] uses Markov models to a great extent in his analysis 
of various social systems. Thonstad [1968] developed a 
perky type model for use with the Norwegian educational 
system. More recently, Branchflower [1970] utilized a Markov 
type model to examine the promotion and appointment policies 
of a university and the resulting rank distributions of 
faculty. Branchflower assumed the promotion policies were 
fixed and used his model to predict the faculty distribution 
in the next time period using the current faculty distributions 
and actual appointments to each rank. In his analysis, 
Branchflower defines three sets of rank distributions: (1) an 
Attainable set, (2) the Maintainable set, and (3) the 
Containment set. An Attainable set is a set of rank distri- 
butions which can be reached from some previously attainable 
set. The Maintainable set is better known as the Steady 
State set of rank distributions. The Containment set is the 
set of rank distributions which.can be returned to infinitely 
often by means of a sequence of appointment policies. 

The analysis in Branchflower [1970] leads us to pose 
several questions about the rank distributions given some 
promotion policy: 

ik Can the present rank distribution be maintained for 


all time by means of appointment policies alone? 


De Can another given rank distribution be attained from 
the present distribution in a finite time period, or can the 
desired distribution never be reached? 

ore What are the rank distributions which can be returned 
to at some time? 

In this paper, our purpose is first, to define explicitly 
the sets of rank distributions of interest associated with 
the model, and to demonstrate various properties of these 
sets. Second, we shall use these definitions and properties 
to develop test procedures to determine if a given distribu- 
tion is an element of a particular set. This analysis will 
provide answers to questions 1 and 2, and will permit us to 


give crude bounds as an approximation to answer question 3. 
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II. THE BASIC MODEL AND NOTATION 

Consider a discrete time Markov model with n+l states 
where state j is the rank of a member of an hierarchical 
Organization. The set of states can be partitioned into two 
subsets: the first n states are transient states and are 
Called active ranks; the second subset is the n+1St state 
which is an absorbing state. Let Pi be the fraction in 
state j at time t+l that were in state i at time t. Then 


n+l 


x op 
jel 


Ae =o a=) 2 ae ee, - 


We represent these fractions in a transition matrix 


P W 
QO = 
0 1 
where P is ann xn matrix, [pis]. tp JH ly 2yp ees Dy and. Was 


an n dimensional column vector called the loss vector with 


W: P is called the Promotion matrix and describes 


fee §i n+1° 
the transition process of the active ranks in the organization. 


n 
It has the property that 12 
j=l 


Pi3 = a i=], 25... en Tar lt la SErlot 


inequality holding for at least one row. 

We assume that the organization is in an equilibrium 
condition, and thus keep the total number in the organization 
constant. We normalize by assuming this number is one. 

Let x(t) be an n dimensional row vector, x(t) = (x, (t), 


Xo (t),..-,Xp, (t)) where x. (t) is the fraction in the ith rank 


idl 


n 
at time t. Note that 2% x; (t) ee el 0 -e bet fe eeebe 
i=l 


an n dimensional row OCC f(t) = (Skt) « caeteets where 
£;(t) is the fractional number of new appointments to rank 
iat time t, ieee) > 0. 

The rank distribution at time t+l is dependent on the 
rank distribution at time t, the Promotion matrix, and the 


appointment policy at time t+l. In vector notation, 
x(t+1) = x(t)°-P + £(t+1). (1) 


The requirement that the total number in the organization be 
constant requires that the total fraction appointed to the 
active ranks be equal to the total fraction lost to the | 
an e be 
absorbing state: 2 £, (t+1) = x(t) ew. (2) 
i=1 

We use the following notation throughout the remainder 

of the paper: eis ann dimensional vector with all components 


equal to 1; e. is an n dimensional unit vector with the jth 


j 
component equal to l. It is apparent by the context whether 


they are row or column vectors. 
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Tit. SET. DEFEND Tons 


The set of all rank distributions, F, is defined as 
Bes! (oc; cies ce 0 eee 


Note that the time notation is omitted when it does not affect 
the definition. 

An Attainable set is a set of rank distributions which 
can be reached in one step from a previous set of rank distri- 
butions. Let m be the number of steps (time periods) from 


the set of all rank distributions, F, then 


2 WS oe pe AUS Ses IME LAT Ie A-1) PaNinisi= sD a eal 


Consider the sequence of Attainable sets, {A,}. If the 
sequence converges to a set C we call the limiting set the 
Containment set. Mathematically, we have 


C = lim 


mse @ 


bal Gl bee oe CAM AUI2 mod (3 19 485 ae 


msaseo 


LekR pes Aci ee a 4 ee) Mabe cle) Mabie Graleh see 


This definition is somewhat unusual in that the Containment 
wees defined in terms of 1tselfr? 

The set of all distributions which can be maintained in 
succesSive time periods is called the Steady State set, M. 


Let f(t) = f for every t, then since x(t+l) = x(t), from 


3 


(1) we have x(t) = x(t)°P + £ 


x(t)°(I-P) = Engei0 
so we define 
Me=miocec (1 Pee cee Fy ee 


If a rank distribution can be reached from the Steady 
State set, M, then we call a set of such distributions a 
Reachable set. All distributions which are reachable from M 


Paem steps are in the set B,, where 


= le ee) Ze Bea! , m= 1,2,... and 


Bo = M. 


Now consider the sequence of Reachable sets, {Bt}. ilee 
this sequence converges to a set B, we call the limiting set 
the Reached set. Mathematically, we have 

Be= | dame, 

mws°2o 


AEE Tigo eae 2 Pee Cte eZee eB 
maw7aso 


la, ence ete ez ce 9b) a 2 teethe Slane xi sits. 


I 


Thus by their definitions, we see that the Containment set and 
the Reached set must be the same set when they VeOchre xX1sS tame aia 


section IV we prove they always exist. 
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IV. SET PROPERIEES 


In this section, we list several properties of the sets 


defined in section III without proofs. We then prove several 
lemmas and two Theorems which establish the existence of the 
Containment set and Reached set. 

PROPERTY 1: The set of all rank dis tir LOWE Onis, al eas 


convex and closed. 


PROPERTY 2: Any Attainable set, ee 1S convex and closed. 

PROPERTY 3: The Steady State ace M, 1s convex and 
closed. 

PROPERTY 4: Any Reachable =o Ba! 1s convex and closed. 


LEMMA 1: An Attainable set, An? is a subset of every 
previously attainable set. 
Proof: F = Ag = {x; xre = 1, x > 0} 
A, = tii Xx 2 zie Poe nemeaz: iE Ag} 


. = = 
Pick x € Ay x <6 Ao therefore Ay = Ag = F 


Assume A,_j S 1, = 


eee = Ay & Ag- 


Define A, (2) Sx 2 Zee ke) eee eso 
A = - ; A(z) + Pick z* € A,_j, then since A,_, $ An-2' 
m-1 


LRG oho PAIS 5 ID. Pi = A,-2! rani 


A,, (2*) = An-1° 


-_ 

N 
oe 

— 
iI 


So An = U A, (2) = ALi: 
ZEA 
m~-1 
LEMMA 2: The Steady State set, M, is a subset of any 


Attainable set, An’ 


1S 


Proofs: M = (x7 eRe oc FH}. MS Ay =F byedek inageiom, 
= 1 * = * z > keep, 
Assume M A-1° PAC he Zeer z* € ALy andeZ 2 aaa 
Therefore, z* € A, by definition of A,, and M* Ap. 


LEMMA 3: Any Reachable set, B,, contains every previously 
reachable set. 
ELOOL at! —lo nme fx > Ee, xX € F}. Trivially, By = By: 
Assume that Bpn-2 = By) where Bm-1 = {x; x > z°P, 
ape Ppa ce Eee 
Define the following sets: B,_,(z) = {x; x > a2-P, 


ANGE CZ)ee= ane xe moe F}', 


gic Calida Zu Ci ae, 


Zz & Bp-)- Piek a ¢ Bee Since Bes, | B81 ,02e eee 


SO) ende Bey (2) — Bata? eo inlces2 2S vae Berane 
SENG (By) ee om auras Br 1 DY eoevens (ea) UB, (2) 
zZ€ By-y ZE Ba-l 
= B.. 
m 


LEMMA 4: If w: > O for at least one i and n > 1, then 


at 
the Steady State set, M, contains an infinite number of points. 
EOo bf? 
(i) Non-empty. Assume that M is empty. Then there 
Ges Otc cieonee SUCH Behotewe (iF 0. a Bocause 
n 


) 
j=1 


Pi; <1, i=1,2,...,n with strict inequality holding 


for at least one 1, (r-p)~1 exists with all elements 
non-negative. Let N = (I-p)-1. «If M is empty, then 


there does not exist an x such that x = £°N > 0. But 


*n.. and since f. nO wand Ne On fOr sey ozs, 


n 
Xs = Lb £f: 
i= 1) 


al 
J il 1) 


ae x; > 0 and therefore M is non-empty. 
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(ii) Infinite Interior. Assume that M contains only 
one point, then x(t + 1) = x(t) for every t and any 
appointment vector. We rewrite equation (1) as x(t) = 
x(t)*P + £(t + 1) = £(t + 1)-(1-P)~t. We also require 
that f(t + 1)-e = x(t)ew. We can construct the 
appointment vector in the following manner: let 

fat al) 


+ Te oa peCumevne Ven CS) x (Aan 


= one | O55 an ad 


Then f(t + 1) = (04 1%5 720070.) and 
(02) SARE jo oo e heme) (3) 


Consider the appointment vector f(t + 1) when 

a, = x(t)*w and a5 = 03 = ..- = a = 0, then x(t) equals 
x(t)*w multiplying the first row of (r-p) 7}, Since 

any f which solves (3) generates a point in M, consider 
the appointment vector f(t + 1) when a5 = x(t)*w and 

Qp = A, = -2- = AQ = 0, then x(t) equals x(t)*w multiply- 
ing the second row of Ge) ack Since x(t) is one point, 
and x(t)*w is a scalar, the first two rows of (r-p)~1 
are identical. This implies that r((I-p) 7+) < n and 
that ((r-p) 1) 71 does not exist which is a contradiction. 
Thus since M is non-empty, M contains more than one 
point. Since M is convex and closed, it contains an 
infinite number of points. 


THEOREM 1: The Containment Set, C, and the Reached Set, 


B, exist, are identical, and non-empty. 


iy 








Proof: Consider the sequence of sets {A}. By Lemmas 
1 and 2, M= ...5 A, = Se cleo SA, SF. By Lemma 
4M +: o> i = $6, n=1,2,.-. .- Thus we have a mono- 
tone decreasing sequence of sets which is bounded below, 
therefore the limiting set soe exist and iS non-empty. 
Now consider the sequence of sets ibe where Bo = M. 


B GS... SF. 


By Lemma 3, M © BgS BS .-. SBS . 


Thus we have a monotone increasing sequence of sets 
bounded above by F. Therefore a limiting set exists. 
As noted in section III the Containment set, C, and 
the Reached set, B, are the same set. 

THEOREM 2: The Steady State set, M, any Attainable set, 
as the Containment set, C, any Reachable set, Bre and the 
set of all rank distributions, F, are identical if and only 
if the Promotion matrix, P, is diagonal. 


Proof: Since M & Ee Cc SA, SCF, we need only show that 


Me= F. 

(1) Assume that P is diagonal. We can write the 

distribution vector as x(t) =a e, tae +...tae, 
ae 22 n 


where Os Be) see lbse BO; 1 then x(t) ¢ F. Rewrite 
equation (1) x(t + 1) = x(t)°P + £(t + 1). We require 
f(t + 1)-e = x(t)*(I-P)-e. Since P is diagonal, x(t)° 
(I-P) is a suitable choice for f(t + 1), (£, (t ag dl) eae), 
So we can now write x(t + 1) = x(t)*P + x(t)*(I-P) = x(t). 


Now x(t + 1) = x(t) and x(t) ¢€ F>? x(t) € M, . oe. Fb Sarre 


And since M& F, M=F. 


ire 








(i1) Assume M= F. x(t + 1) = x(t) for all t. 
Assume that P 1s not diagonal. x(t) = x(t)*P + £(t + 1). 


Since x(t) ¢ F, we can write x(t) = + + ... + 


“sie.” Same 
nen, A; > Ou ria = 1. In particular, let oo lig 


a, = 0, i¢ 1. So x(t) = a ,e, and (1,0,...,0) = 


(1,0,..-,0)* | Py, Pip +++ Puy 


° ee ca 


Phi a ses) Pon 


= (Pyqr Pyo rece Py) + (Ey fo ree, f,): This requires 


that the following equations be satisfied: 


+ = 
eT) Pil : 
f. + Pia = 0 
ry ile ; 


J 


repeat this procedure for a; = 1, a, = 0, k eS 


2,3,...,n. Thus py3 = 0, i + j and P is diagonal. 


Since go 2 Ge ee 0, we have f. = Pi = 0, j +1. We 


1h 





V. EXTREME POINTS OF THE VAREOGUS SETS 

The characterizations of the various sets in the model 
in the above section do not describe them in a useful manner. 
If we knew the extreme points of all of the Sets, then we 
would have almost all of the information about the sets, 

Since any interior point of a set is simply the convex combin- 
ation of its extreme points. 

For two sets in particular, the extreme points are easily 
obtained. The set of all rank distributions, F, is the 
manaamnental samplex in E, with extreme points, e;, i= 1,2,...n. 
The subset of F which is most easily described is the Steady 
State set, M. Each vector in the set satisfies the inequality 
ie ee P. This can be written in equation form as x*(I-P) =f 


where f is the vector of appointments to the different ranks. 


We define N = (T-p) 71 so we write x = f*N. Since we require 

Peemee lt and Since £ can be written £ = Ya; (x*(I=P)-e), G- 0, 
i 

ta. = 1, we write the equation in an equivalent form as x = f£*B 


i 





New 
where f = Lae; and Bi5* 5 -J , That is, the extreme points of 
i eshig 
a 


the Steady State set are given by the rows of (r-p) “1 after the 
rows are normalized. 

The problem of determining the extreme points of the other 
sets in the model is much more difficult. If we are working in 
three dimensions, it is not too difficult to calculate all of 


the extreme points of one of the sets, A, or BW! for a small m. 


i 
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Each set in the model has the property that each point in the 
set generates a subsimplex in the fundamental simplex. The 
extreme points of the generated set are found from the convex 
hull of the subsimplices generated by the extreme points of 
the original set. In three dimensions we use a geometrical 
method to determine the extreme points of the convex hull of 
the generated subsimplices. If we know the extreme points of 
Peeeor some m, then we project A, onto a plane generated by 
two unit vectors. Then we calculate the extreme points of 
the subsimplex generated by each extreme point of A, and project 
these onto the plane. Then we simply test linear combinations 
of the generated extreme points to determine that there are 
no points which lie beyond the linear combination in each of 
the unit directions. The points which satisfy this criterion 
are the extreme points of the set Ay4y). A simil#r procedure is 
used for B,. 

When P is upper triangular, some of the extreme points of 
A, are also extreme points of Any , (e.g. €, is an extreme 
pemme Of A, for every m.) These remain extreme points for all 
mueeceding sets generated from Aj,,]. Similarly, B,, has 
extreme points common to all B, which are extreme points of 
A, as well, and therefore are extreme points of the Containment 
set, C. Figure 1 depicts several attainable sets drawn on 
barycentric coordinates for a three dimensional rank 


distribution where 
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FIGURE 1. ATTAINABLE SETS ON BARYCENTRIC COORDINATES 
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Table I contains the coordinates of the numbered points on 
Figure 1. The heavy line in the figure represents the known 
portion of the boundary of the Containment set after calculating 
A7. The portion of the boundary between points 9 and 10 

Changed very little between Ag and A> so the boundary of Av 

is probably a good approximation of the boundary of the 


Containment set in this case. 


TABLE I: POINTS ASSOCIATED WITH FIGURE 1. 


; dh i 
“a 

39, 73/7231 2.35269 | 226615 176445! -0882225 noe 
ABR 229401. 5259321-49206 1.500700) 5355714 


BSED .3337] .32179].32179|.322849 -3762061| .5]1/3 | 


It is not clear how this graphical projection method can 
be extended to higher dimensions to obtain the extreme points. 
Nor is it clear, based on the above example, whether or not 


the Containment set has a finite number of extreme points. 


Z3 








VI. THE DEVELOPMENT OF SEVERAL Tes Es 

Since it is not feasible to analytically determine the 
extreme points of all of the sets, we must devise some method 
to answer our original questions. The answers to those 
questions depend on knowing the extreme points of all of the 
sets. To provide partial answers to the questions, we develop 
test procedures to determine if a given vector is an element 
of a particular set, and to determine if a given vector can 
be returned to inm steps or reached from another specified 
vector inm steps. We know that a given vector cannot be 
returned to at any time if there exists an i such that the 
feot vector is not in A;. 

Fach transition from one set to another reguires that a 
certain set of inequalities be satisfied. When we consider 
the system of inequalities that account for all of the transi- 
tions, we need only test for feasibility. If no feasible 
solution exists, then the test vector fails the test. Several 


specific tests are developed in the next sections. 


A. Ay bio FOR Bn 
In this section, we develop a test to determine if a given 

rank distribution is an element of the set eS enecatl that 

Mee 12 5 2, > Znny Pe 2 


2B, 43) (SE 7 dgtelo jel 2s ae 


m m 


and that B,, is the set of rank distributions which is reachable 
in m steps from the Steady State set. Row vectors have been 
used in all developments. To simplify notation, the test 


inequalities will be specified in tableau form and vectors 
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will be column vectors. Transposition is indicated for the 
promotion matrix where it appears. 
We now write the sequence of inequalities which must be 


Peielomled if thesyectom xX 15 limene Set B. as 


0 . 4 Zy 2 2s Z Z 


a5, 
ac 
l 
be 
l 
= 
[A [A 
ro) S 


Ae) 
-j 
j 
IH 
[A 
ro) 


yg 
| 

LA 

~ 


fj 2 oo fH 


e = |] 


The constraint matrix appears deceptively smaller than it 


is. Since each z; is ann vector, the dimensions of the matrix 


1 
are mn + m + n rows and mn columns. Even though the constraint 
matrix iS very sparse, we were unable to apply any algorithm 


which would process the problem more efficiently than the 


Simplex algorithm. 
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It is clear that a feasible solution to (4) 1S simply a 
Phase I Simplex solution. The method which we use to obtain 
ta2s solution 1S to add an artificial variable toveach or ene 
feasibility constraints, the last m rows. We then minimize 
the sum of the artificial variables. A feasible solution to 
the new system, (4) plus the artificial variables, always 
exists. The result of the test is determined by the values of 
the artificial variables. If all of the artificial variables 
are equal to zero, then (4) has a feasible solution and we 
eeuclude that the test vector, X, 1s an clement Of BL. If 
any of the artificial variables are not zero, then a feasible 
solution to (4) does not exist and we conclude that X is not 


in Bae 


B. A TEST FOR A, 
In this section, we develop a test to determine if a 
feeven rank distribution is an element of°A,, 1-e. 2f it can 
be reached in m steps from some point in F. Recall that 
eee 12; 7A NA ae Dyes Sue A,-]} for m=1,2,... 


We write the set of inequalities as 
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pt < xT (5) 


FI ooo 


The constraint matrix has mn + m rows and mn columns. We 
now test (5) for a feasible solution in the same manner as 
miemeceot fOr BW. The results of the test are anterpreted in 
fiemoame manner as those for B,. Note that the constraint 
matrix in (5) is the same as that in (4) without the first n 


rows. 


eC A TEST FOR RETURN IN m STEPS 

In this section, we develop a test to determine whether 
or not a given rank distribution, X, can be returned to inm 
steps. If X is an element of the Steady State set, we know 
that it can be returned to in one step, -i.e. it can be 


Maintained forever. So we are interested ina test for m 
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greater than one. Since there are m transitions, there must 
be m-l intermediate vectors because the initial and terminal 
vectors are both X. The first intermediate aS. 24) must 
satisfy Z, > XP. The second intermediate vente Zor must 

Satisfy 29 > Z1°P. This process continues until the last 


intermediate vector, Z,_,, must Satisfy z “=P <eee We Wate 


Go 


the sequence of inequalities in tableau form as 


Z4 Z. - 2 2 Zo Z 4-1 
I > pi.xT 
pt -I < 0 
p* < 0 
oe =I < 0 
(6) 
ptr < xt 
e = ] 


0) 
Il 
eee 


Again we notice that there is a great Similarity between 


Pieroni Om thls tableau and  tnoSe ron the tests LOK he and 
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B The first n inequalities in (6) are different from those 


eye 
in (4), as well as there being one less intermediate vector. 
The constraint matrix in (6) has mn + m +n rows and mn - n 
columns. 

As with the tests for A, and Bae any method of finding 
a feasible solution to (6) is acceptable. One method would 
be to add artificial variables not only to the last m rows 
of (6), but also to the first n rows, then minimizing the 
sum of the artificial variables. If all of the artificial 
variables are zero, then we conclude that the rank distri- 
bution, X, can be returned to inm steps. The method used 
in section VII involved adding artificial variables only to 
the last m rows of (6), then minimizing the sum of the artifi- 
cial variables. This system may not always have a feasible 
solution. The test requires that if the system has a feasible 
solution, then if all of the artificial variables are zero, 


we conclude that X can be returned to inom steps. 


Dy A TEST FOR CHANGING DISTRIBUTIONS IN m STEPS 

We now develop a test to determine, given some distribu- 
tion vector, Xi whether or not we can reach another specified 
distribution vector, Xz, in m steps. This is nearly identical 
to the problem discussed in the above section, i.e. we have a 
fixed initial vector and a fixed terminal vector. The con- 


straint matrix is the same as (6) but now the right hand side is 
T 
a! 


the same as the test for return inm steps. 


ee Ore pk 1,---, 1]. The test procedure is then 


Has) 





Le A TEST FOR EXTREME POINTS OF THE, CONTAINMENT SET 
The similarity between the tests for A, and B, noted in 
sections A and B above permits us to develop a procedure with 


which we can determine extreme points of the Containment set, 


om 
Let G) = Hare), O,---, 0] be ann by mn matrix and let 
pr =I 
eo aot 8 
pr -I 
pi 
= 
e 
-e 
e 
-e 
e 
be an mn + 2m by mn matrix. Now let Xo = [One . .005e): eee 
be an mn + 2m by 1 column vector and Z@ = ee i be an 
m- 


mn by 1 column vector. Then the test for ee reduces to finding 
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a feasible Solucvon to 


‘2 < X - (7) 


G 
2 


Zz 


If there is a Z* which satisfies (7) then the test for By 
would seem to result in testing to aoe if G,°4* 20. = HOWeW eta, 
there is no guarantee that Z* is unique and there could 

exist a Z' which solves (7) and also solves G,°Z' > 0 for the 
same test vector X. So this would be a) untied eter test For 
inclusion in B,, but not necessary. The only necessary and 
sufficient test is that described in section III.A. 


Now let H, = [(I-PT),0,...,0] be an n by mn + n matrix 


and let 
pr -f 
pes 
; =F 
. pr -x 
= 
= 
Be = 
= 
—— 
e 
=e 


=e 


cal 





be an mn + n + 2m by mn + n matrix. Let Z‘' = 125 ae be 

an mn + n by 1 column vector and let X4 = (ieee oe 

be an mn + n + 2m by 1 column vector. This amounts to letting 
the test vector, xt become an unknown, Ze Then all of the 
basic feasible solutions to H5°Z' < are extreme points of 

A: If a point Po eal econo Pel eroten sana Ba Simultaneously, 
then it is an extreme point of the Containment set, C. Thus 

we need only test each basic feasible solution to see if 

Mee > 0. IL it is, then 4’ is san extreme point of 


Due to the computational difficulties associated with 


this test, it is presented for theoretical interest only. 
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VII. <A CASE STUDY 

The tests developed in section VI above were applied to 
faculty distributions in the College of Engineering at the 
University of California at Berkeley based on data from 1960 
to 1968. The bulk of the data is in Branchflower [1970] and 
only that which is relevant is reproduced here. 

We used the Mathematical Programming System on the IBM 
360/67 computer at the W. R. CHURCH Computer Center, Naval 
Postgraduate School to conduct most of the tests. This is 
really a trivial exercise once the data cards are prepared 
for the control programs. One FORTRAN IV program was written 
to generate the data cards for the test for B,. The same 
eards were ised sto test tor A sab tem mcmoving lec aid omc ise 
taining to the first n rows of the constraint matrix. Another 
FORTRAN IV program was written which was used to generate 
data cards for testing return inm steps. Again, hand modifi- 
cation of the right hand side permitted the same cards to be 
used to test for changing distributions inm steps. 

One other FORTRAN IV program was written which was used 
to determine the extreme points of the steady state set. It 
was also used to test to determine if a given vector was an 
element of the Steady State set. 

All of the tests were verified for accuracy by using a 
test case, namely the three state promotion matrix in section 


V and points which were in known sets from Figure l. 
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FIGURE 2. PROMOTION MATRIX FROM 1960-1968 DATA 


TABLE II. TEST VECTORS 


<j) (EE 5)5))) 2 GALS SS | X3 (1967) X4 (1968) 


-0161427 - OS 4639 a ~004/ 612 
-US02Zi32 -0567010 nU62G 012 .0190476 
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-0374331 e025 7751. .0628019 .0761904 
70934759 -VGLva56 S(O ibis SiO) 6 sUaZzo0Ug 
-0481283 -0463917 ONS 7215 OE) . 0666666 
-1390374 s097 Wisc -UiGza 2510 -0952380 
S64 7.0 -1494845 - 1400966 -1047619 
ot 209'4 s11340;20 seo 420 siZ230 005 
- 1390374 710309 27 -1062801 -0904761 
,UI6Z 5016 Pelgilycyayey oy) ~1014492 O47 Gers 
- 0655614 .0670103 .0676328 -1047619 
20213903 -0463917 Foye face .0761904 
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For the case study, we used the thirteen active states 
of the faculty distribution and the 1960-1968 promotion 
matrix from Branchflower [1970]. This is reproduced as 
megure: 2. The actual distributions of faculty, formige> oc 
were normalized and used as the four test vectors. These 
are included in Table II. 

The test vectors were tested for being elements of ‘Ge 
and Bar mM = 1,2,...,/7, and were tested to determine if they 
could be returned to inm steps for m = 2,3,...,8. The 
resulting location of the vectors of rank distributions on 


the fundamental simplex is indicated in Table III. 


PASE Tif. TEST RESULTS 





Recall that if X; E A. then after one step, X > X.°P, 


i+] 
X.4, 1S at least an element of Assi: The results in Table III 
are inconsistent with this theoretical result and require an 
explanation. There are two reasons why this could have occured. 
The first reason is that the Promotion matrix was an average 
over an eight year period and did not exactly represent the 
true Promotion matrix for each year. The second reason is that 


the faculty size did not remain constant, but in fact was 


increasing. An expanding organization can appoint personnel 
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to replace the losses and then hire the additional number 
for expansion into any state. When P is upper triangular, 
the extreme points of successive a may tend to move away 
from the lower ranks. For a three state asRenTD Ue, consider 


an initial rank distribution of (1, 0, 0) with 


.S, 24 0 
Pp = 0 ALS aS: 
0 0 <7 


then after one transition, the rank distribution before new 
appointments is (.5, .4, 0). Say the number of persons in 
the organization is constant at 100, then the distribution of 
persons in the ranks is (50, 40, 0) with a loss of 10 to 
replace. Appointments could result in the distributions 
(60, 40, 0), (50, 50, 0), or (50, 40, 10) or any linear 
combination of them. If the organization expands by 10 persons, 
then the distributions could be (70, 40, 0), (50, 60, 0), or 
(50, 40, 20) or any linear combination of ti.em. If we 
normalize, then these expanded distributions are outside A, 
for the constant size case. The result is that the possible 
distributions for the expanding size model can be outside of 
the normal set of distributions obtained for the constant size 
model. 

We then aggregated several ranks in the model to determine 
if the aggregated model behaved in the same way. The first 


four ranks are Assistant Professors, the next three ranks are 
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Associate Professors, and the last six ranks are Full 


Professors. The resulting Promotion matrix is 


.7649 .1752 0 
Pp = 0 .8030 .1841 
0 0 .9707 


The normalized rank distributions for 1965-1968 are 


listed in Table IV. 


TABLE IV. TEST VECTORS: AGGREGATED MODEL 


X, (1965) X, (1966) X, (1967) X, (1968) 


2S GG -195876 207022 ee 952 
- 240642 - 206186 IS 237 ae Lai2eo 


Se, he e271 936 779034 -604762 





These vectors were tested to determine their location on 


the fundamental simplex. The results are shown in Table V. 


TABLE V. TEST RESULTS: AGGREGATED MODEL 





a7 





These sesults Using sehevaggqreqateds mode) warcurs a1) iia e 
the results of the thirteen state model. We noted above that 
as faculty size increases, the rank distribution may move to 
a superset of the set containing the original distribution. 

In fact, this happens in 1966 for ain GmoRaewe of seven faculty 
members (2072), and also in 1967 for an increase of thirteen 
faculty members. In 1968, there was an increase of three 
faculty members, but the rank distribution moved to a subset 


of the set containing the 1967 rank distribution. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 

We have established several properties (of ssevereceds 
of rank distribution in an hierarchical organization which 
have not been previously examined. In Seer, based on 
these properties, we developed several tests to determine if 
a specified rank distribution is an element of certain sets. 
The tests were applied to data from a case study of the 
College of Engineering faculty at the University of California, 
Berkeley. 

Several questions remain unanswered and are possible 
areas for further study. Is there some neighborhood of a 
given rank distribution which can be maintained? Is there 
an efficient algorithm which will solve the test inequalities 
other than the Simplex method? Is there a method of reducing 
the number of basic feasible solutions which must be found to 
determine the extreme points of the Containment set? Does 


the Containment set have a finite number of extreme points? 


be 
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